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ABSTRACT 

Radiative transfer in planetary atmospheres is usually treated in the static 
limit, i.e., neglecting atmospheric motions. We argue that hot Jupiter atmo- 
spheres, with possibly fast (sonic) wind speeds, may require a more strongly 
coupled treatment, formally in the regime of radiation- hydrodynamics. To low- 
est order inv/c, relativistic Doppler shifts distort line profiles along optical paths 
with finite wind velocity gradients. This leads to flow-dependent deviations in the 
effective emission and absorption properties of the atmospheric medium. Eval- 
uating the overall impact of these distortions on the radiative structure of a 
dynamic atmosphere is non-trivial. We present transmissivity and systematic 
equivalent width excess calculations which suggest possibly important conse- 
quences for radiation transport in hot Jupiter atmospheres. If winds are fast 
and bulk Doppler shifts are indeed important for the global radiative balance, 
accurate modeling and reliable data interpretation for hot Jupiter atmospheres 
may prove challenging: it would involve anisotropic and dynamic radiative trans- 
fer in a coupled radiation-hydrodynamical flow. On the bright side, it would also 
imply that the emergent properties of hot Jupiter atmospheres are more direct 
tracers of their atmospheric flows than is the case for Solar System planets. 
Radiation-hydrodynamics may also influence radiative transfer in other classes 
of hot exoplanetary atmospheres with fast winds. 

Subject headings: planetary systems, radiative transfer 



1. Introduction 



In recent years, it has become possible to remotely observe the atmospheres of some 
exoplanets found in close-in orbits around nearby stars. Many of these observations, which 
include eclipse, phase curve, and transit photometric or spectros copic measurements, have 



focused on the specific class of exoplanets known as hot Jupiters (|Char 
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20071 ). These gaseous giant planets have high atmospheric temperatures, from the close prox- 
imity to their parent stars, slow spin rotation rates, as inferred from their expected state of 
tidal synchronization, and they are thus subject to an u nusual steady pattern of hemispheric 



insolation, w ith permanent day and night sides (e.g., ISeager et al. 
2008bl . l2010[ ). 



2005; Showman et al. 



A variety of diagnostics about physical conditions in the atmospheres of a few specific hot 
Jupiters have be en presented in the literature, from the presence of high altitude haze (e.g., 
Pont et al.ll2008[) to the existence of ver tical temperature inversions ("stratospheres", e.g.. 



Burrows et al.ll2008l : iFortney et al. 



of radiatively active spec ies (e.g. 



2008), as well as constraints on the chemical abundances 



Charbonneau et al.l l2002t iTinetti et al.l 120071 ; ISing et al. 



20081 : ISwain et al.l l2008bl ). In most cases, these interpretations rely on one-dimensional, 
steady-state radiative transfer models describing in great detail the atmospheric structure 
that is expected at radiative and chemical equilibrium. On the other hand, these models usu- 
ally consider only globally averaged atmospheric properties and they often ignore the role of 
horizontal heat transport by advection on the energetic balance of the modeled atmospheres. 
Such shortcomings of radiative transfer models could eventually limit our ability to reliably 
infer the physical conditions present in remotely observed exoplanetary atmospheres. 

Planetary atmospheres which satisfy global radiative equilibrium are generally not in 
local radiative equilibrium, that is along a specific vertical column, because of finite contribu- 
tions from horizontal heat fluxes. The role of atmospheric circulation in shaping the structure 
and properties of hot Jupiter atmospheres h as been recognized as an increasingly impor- 



tant ingredient over the last few years (e.g., iFortney et al.l l2006l ) and several groups have 
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Even in a three-dimensional, time-dependent atmospheric model of the type known as 
"General Circulation Model" (GCM), radiative transfer is traditionally treated in the static 
approximation. In this limit, the coupling between the atmospheric flow and the radiative 
heat transport enters only via the energy equation for the moving atmospheric fluid, so that 
this coupling is purely thermodynamic in nature. The radiative transfer component of the 
problem, which is needed to continuously evaluate the net diabatic heating or cooling rate 
of the gas in motion, is solved independently of the atmospheric motions themselves, as if 
the atmosphere were actually static. 
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Here, we argue that this level of thermodynamic coupling, which is a standard approxi- 
mation for Solar System planetary atmospheres, may be insufficient to describe accurately the 
nature of energy transport in hot Jupiter atmospheres. A more strongly coupled treatment 
known as radiation-hydrodynamics, in which flow velocities enter the radiation transport 
problem explicitly, may be required to describe radiative transfer in hot Jupiter atmospheres 
with fast, sonic or transonic wind speeds. 

The remainder of this paper is organized as follows. In §21 we outline the general 
formalism for radiation transport in a dynamic atmosphere. We motivate the relevance of 
this radiation-hydrodynamics regime for hot Jupiter atmospheres with fast winds in ||3l We 
discuss the potential magnitude of deviations from the static limit for the global radiative 
balance of hot Jupiter atmospheres in §1] and we conclude in §3 



2. Radiative Transfer in Dynamic Atmospheres 



Radiation-hydrodynamics requires a fully relativistic treatrnent to properl y account for 
Doppler shifts and aberration effects (IMihalas fc MihalasI Il984j : ICastorl 120041 ). The time- 
dependent radiation transport equation in the fluid comoving frame, to leadi i ig ord er in the 
ratio of the fluid velo city to the speed of light, f/c, was derived by iBuchlerl ( Il983l . see also 
Hsieh &: Spiegellll976l ). The breakdown of terms in this rather general equation, which may or 
may not be discarded depending on the nature of the physical problem at hand, is a delicate 
matter (jCastoiil2004j : iMihalas fc MihalasI 1 19841 ). In the limit of slow fluid motions {v/c <^ 1) 
and thus fast light-transit times, however, there is general agreement that the leading order 
effect of fluid motions i s simply to induce Doppler shifts in the i nstantaneous radiation fleld 
seen by the fluid (e.g., iRvbickil EgTol : iRvbicki fc Hummeil EgTsl l 1^ This results in modifled 
absorption and emission coefficients for radiation-matter interaction, particularly when these 
interactions are dominated by line opacities. In the remainder of this work, we focus our 
discussion on line opacities, which are expected to dominate in hot Jupiter atmospheres, and 
neglect continuum sources of opacities. 

Much of the astrophysical literature on radiation-hydrodynamics is focused on the ra- 
diative transfer o f radially expanding s tellar winds. We start our discussion closely following 
the formalism of iRybicki fc Hummen ( 119781 ). To lowest order in v/c, the radiative transfer 
equation in a moving medium can be written 



^Effects due to relativistic aberration are thought to be significant only at higher order m v/c ( Castor 
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n-V/(r,n, z/) = -fctotW$ 



{i> - z/q) + — n-w(r) 
c 



[/(r,n, z/) - 5(r,n, z/)], 



'11 



where J(r, n, v) is the monochromatic intensity (per unit sohd angle) at frequency u and 
spatial location r in the direction defined by the unit vector n, 5* is the source function, 
v{r) is the material velocity field, fctot is the line integrated opacity, z/q is the line central 
frequency and $(z^) is the line profile function, normalized to unity. The standard radiative 
transfer equation in the static limit is recovered by forcing v{r) to zero. 

Equation ^ expresses the fact that any projected velocity gradient along an optical 
path in the direction n leads to a shift in the line profile, and thus a change in the effective 
absorption coefficient of the medium, in proportion to the Doppler term n-v{r)vQ/ c entering 
the line function $. By contrast, a uniform velocity field leads to a more trivial, uniform 
frequency shift. In the context of stellar wind theory, equation ([T]) is typically solved using 
the approximate "Sobolev method," in the lim it of strongly supersonic gradient flow ra otions, 



for radially expanding shells of material (e.g., iRybicki fc Hummer! Il 9781 : ICastorl 120041 ). This 



is not the limit of interest for hot Jupiter atmospheres. 

The literature on radiati ve transfer in static a t mospheres is vast and the rea der is 
directed to the monographs by iGoody fc Yung] (Il989[ ) , iThomas fc StamnesI (|2002[ ) and iLiou 
(I2OO2I ) for detailed accounts on this subject. To make the connection with that work, we now 
introduce the plane-parallel atmosphere approximation, using z as the vertical coordinate 
in the atmosphere. In addition, we introduce the zenith angle 9 (away from the vertical 
axis, pointing upward) and the azimuthal angle (in the plane normal to the vertical). It 
is convenient to introduce the quantity = cos^, which ranges from +1 going up to — 1 
going down in the plane-parallel atmosphere. We further assume that local thermodynamic 
equilibrium is satisfied and we neglect scattering, as is generally justified for the thermal 
portion of the atmospheric radiation spectrumo The source function S then reduces to 
the Planck function, By{T{z)) = By{z). In that case, the formal solution to the radiative 
transfer equation for the monochromatic intensity at level zn m the direction (/i, 0) can be 
written JCoodv fc Y"uiig||l989l : [rhomas fc StamnesI booi lLioull2002l ): 



^It is customary in the treatment of atmospheric radiation to separate the short- wavelength component 
of the spectrum, which is related to beam-like insolation in the presence of scattering (typically in the op- 
tical), from the long-wavelength component, which is related to diffuse thermal emission from the various 
atmospheric layers (typically in the infrared). We follow this convention here, even though the two spec- 
tral components are not as clearly separated for hot Jupiters as they are for "cool" Solar System planets. 
For simplicity, we focus our entire discussion on the long-wavelength thermal component and postpone a 
discussion of possible complications with the short-wavelength portion of the radiation spectrum until 53 
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I„{zo;iJ,,(f)) = I^{0;ij,,(j))Tu{0,Zo;fi,(f))+ / B^{z; /J,, 



dTu{z, zq] fi, I 
dz 



-dz. 



(2) 



where the monochromatic notation has been changed to Ii, and z = corresponds to a 
bounding region for the atmosphere (e.g., the bulk interior), with a one-sided radiative 
boundary condition specified as lu{0; fi, (p). The first term on the RHS of eq. ([2]) describes 
the cumulative absorption of I^{0] fi, (j)) from z = to zq along the optical path, while the 
second term describes the cumulative Planck emission, and its absorption, integrated over 
all layers along the optical path, from z = to zq. The monochromatic transmissivity 



Tu{z,zo;fi,(j)) = exp 



zo 



K{z'; n, (j))p{z'; jj,, (fydz 



(3) 



encapsulates the absorption (and emission) properties of the medium. It is the integral of 
the mass absorption coefficient ky weighted by the amount of absorber (with mass density 
p), between vertical levels z and zq, along the optical path defined by the direction (/x, 0). 

According to equation ([T]), bulk Doppler shifts in a dynamic atmosphere modify the 
radiative transfer by changing the monochromatic absorption coefficient (unit of cross section 
per unit mass). 



•-tot 



{u - uq) + —n-v(r) 

c 



(4) 



or equivalently by changing the monochromatic transmissivities, %,{z, Zq] in equa- 

tion (|2]). If bulk Doppler shifts are important and the material velocity field is anisotropic, 
so is the radiative transfer. 

By property of isotropy of the Planck source function, a large fraction of the thermal 
atmospheric radiation field is carried along rays which are significantly slanted relative to the 
vertical, even though net fiux exchange between atmospheric layers occurs vertically. The 
next step in deriving vertical fiux equations for the atmospheric radiation problem typically 
involve s separate angular integrations for the ascending (m > 0) and descending (/i < 0) 
fiuxes JCoodv fc Yunejll989l : iThomas fc Stamnesl[20o3 : lLioull2002h . of the type 



?± 



2n 



±1 



Iy{jj,, (fyjjLdjJidcl). 



(5) 



To account for the phase space available at slanted angles, it is customary in "two- 
stream" formulations to replace all angular integrals over transmissivities by a single average 
transmissivity value for a characteristic zenith angle, /xq = cos 6*0: 
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2-K /-il 



^0 

TT exp 



exp 



ku{z'] /i, (j))p{z'; jj, (f))dz' 



fidfidcp 



= Tu{z, Zq), 



(6) 



where %, is the hemispherically-averaged "diffusive transmissivity." This approximation has 
been shown to yield errors < 1.5% in ty pical apphcations for static atmospheres, when a 
diffusivity factor 1 //in = 1 .66 is adopted (IRodgers fc Walshawl Il966l : iGoody &: Yung Il989l : 
Thomas fc Stamned 120021 : iLioul 120021 ). which corresponds to a zenith angle 6q = 53 deg. 



This relatively large value of the effective zenith angle illustrates well the large phase space 
available for thermal radiation at slanted angles. It indicates that the projected wind velocity 
along a typical optical path for thermal atmospheric radiation can be a significant fraction of 
the full horizontal wind speedjfl even though net radiative exchanges between the atmospheric 
layers occur in the vertical. 



3. Hot Jupiter Atmospheres 

To help us focus our discussion further, we now turn to issues specific to hot Jupiter 
atmospheres. 



3.1. Projected Wind Velocities 



Over the last few years, various circulation models have indicated that wi nd speeds could 
reach sonic or even supersonic values in the upper atmospheres of hot Jupiters (jCooper fc Showman 



2005, 2006; Dobbs-Dixon fc Lin]l2008: Showman et al.ll2008a. 2009; Rauscher fc Menou 



2009). 



Here, for concreteness, we use the specific hot Jupiter model described in lRauscher fc Menou 
(120091 ) to evaluate the magnitude of projected velocity gradients along representative optical 
paths for thermal radiation in a dynamic hot Jupiter atmosphere. 

A specific optical path is defined by a starting location at the bottom of the model 
atmosphere and by a unit vector n which defines the path orientation. We calculate the 
projected wind velocity along this path throughout the entire model atmosphere as 



■^Large scale motions are predominantly horizontal in an atmosphere. 



- 7- 



Vpro] = n-Vh = rih-Vh sin6'o, (7) 

where Vh is the wind horizontal velocity vector according to the circulation model, rih is the 
horizontal unit vector projected on the sphere (i.e., defining the north/south and east/west 
directions) and the representative zenith angle 60 = 53 deg is uniformly adopted in all our cal- 
culations. A detailed calculation would require us to take into account the three-dimensional 
geometry of the problem, with varying pressure levels in each of the model vertical columns 
crossed by the slanted optical path. Rather than performing delicate three-dimensional in- 
terpolations between various model columns, we use the profile of wind velocities in a single 
vertical column. That is, we use values of Vh as if the optical path were exactly vertical, even 
though the calculation assumes a zenith angle ^0 = 53 deg and various azimuthal orientations 
for the projection. While this approximation clearly emphasizes vertical velocity gradients 
over horizontal ones, it still captures representative changes in Vh, both in magnitude and 
direction, along the selected optical path (with a specific azimuthal orientation). It should 
thus be sufficient to evaluate the typical magnitude of projected velocity gradients along 
representative optical paths in the model atmosphere. 

Figure d] shows profiles of wind velocity projected along representative optical paths, as 
a function of pressure, p (a proxy for height in the atmosphere). All velocities are expressed 
in units of the local value of the adiabatic sound speed, Vproj{p)/cs(p ;T). The same H2- 



dominated atmospheric gas parameters as in iRauscher fc Menoul (120091 ) are used to calculate 



the sound speed: Cg = Vt^T, with an adiabatic index 7 = 1/(1 — k), k = 0.321 and a gas 
constant TZ = 4593 J kg~^ K~^. The various panels show profiles at the model substellar 
point (a), antistellar point (b), west equatorial terminator (c) and north pole (d). In each 
panel, the various curves show projected velocity profiles for optical paths oriented to the 
east (solid line), north-east (dotted), north (dashed) and north-west (dash-dotted) on the 
sphere. Projected velocity profiles for other cardinal directions can be deduced by symmetry. 

Figure [1] reveals significant gradients in projected wind velocity as one crosses the atmo- 
sphere along representative slanted optical paths for thermal radiation. Velocity differentials 
over one pressure scale height easily amount to ~ 0.2- 0.5 Cg and they exceed Cg in some 
cases, especially high up in the atmosphere. Velocity differentials > Cg are typical when 
crossing several pressure scale heights. 

The different panels in Figure [1] illustrate the diverse character of projected velocity 
profiles at various locations around the planet. Furthermore, panels a), b) and c) exem- 
plify the anisotropic nature of these projected velocity profiles. Projected velocity gradients 
are systematically weak in the north (or equivalently south) direction but they can be- 
come significant when the east-west direction is sufficiently sampled by the optical path 
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under consideration. This is easily understood as resulting from the predominantly zonal 
(east- west) nature of wirids in this and other ho t Jupiter atmospheric circulation models 
( iRauscher &: Menoull2009l : IShowman et al.ll2008bl ). It is then clear that the azimuthal angu- 
lar phase space for thermal radiation will be dominated by optical paths sampling significant 
velocity gradients, with unusually low gradient values relevant only for the small fraction 
of all paths that are closely aligned with the north-south direction. As discussed in ^ to 
the extent that projected velocity gradients of the magnitude shown in Fig. [1] impact the 
transport of thermal radiation, the radiative transfer problem will become anisotropic via the 
sampling of a variety of azimuthal and zenith angles, even when the Planck source function 
itself is isotropic (in the fluid frame). 



3.2. Line Shapes and Widths 



Opacity sources in hot Jupiter atmospheres are largely dominated by discrete atomic 
and molec ular lines, particularly in the thermal portion of the atmospheric radiation spec- 
trum (e.g., ISharp &: Burrows! 120071 : iFreedman et al.ll2008l ). Lines are broadened well beyond 
their quantum mechanical width through Doppler shifts associated with the thermal motions 
of atomic and molecular constituents ( "Doppler broadening" ) and through the effect of colli- 
sions of these atoms and molecules with other gas constituents ("pressure broadening"). The 
relative width of Doppler and pressure broadening depends on the local conditions of density 
and temperature in the atmosphere, as well as on the specific radiative constituent under 
consideration. The general shape of a line that is both Doppler- and pressure-broadened, i.e. 
the detailed functiona l form of the line function $ in eg. ([Tl) , is given by the Voigt function 
(ICoodv fc Yi^ll989[ iThomas fc Stamnesl[20o3 : lLioull2002l ). 



Intuitively, one expects the radiation transport in an atmosphere with strongly pressure- 
broadened lines and very subsonic wind speeds to be relatively insensitive to the additional 
Doppler shifts due to bulk flow velocities projected along various optical paths. Indeed, the 
typical thermal width of a Doppler-broadened line is closely related to the sound speed of the 
atmospheric gas (see eq. [10] below), so that bulk line shifts should only minimally deform 
line shapes in the limit of very subsonic wind speeds. This, combined with the dominance 
of pressure broadening, is the main justification behind the use of a static treatmen t for 
radiation transport in Solar System planetary atmospheres (e.g.. iGoody fc YungI 119891 ). To 
evaluate the possible role of bulk Doppler shifts for radiation transport in hot Jupiter at- 
mospheres, it is thus important to evaluate the shapes and widths of radiative lines in these 
atmospheres. 



Various static, globally-averaged radiative transfer studies have established that the ra- 
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diatively forced regions of hot Jupiter atmospheres are found, broadly speaking, above the 



10-1 bar pressure l evel (e.g.. ISeager &: Sasselovl Il998l : ISudarsky et al.l I2OOOI : iBarman et al. 



20051 : liro et al.l 120051 : iSeager et al.l 120051 ). These models also indicate that photospheric lev- 
el^ vary strongly with wavelength in the optical-IR spectral range, as a result of strong 



variations in line opacit ies, from a few bars to < 10 ^ bars typically (e.g., ISeager et al.ll2005 
Sharp fc Burrowsll2007l ). Regions in the pressure range from ~ 10 to 10~^ bars are thus of 



particular interest for the study of radiation transport in dynamic hot Jupiter atmospheres. 

A detailed account of the typical treatment of lines in static rad i ative transfer calcu- 
lations for hot Jupiter atmospheres is provided by ISharp &: Burrowsl (120071 ). The general 
Voigt profile of a pressure- and Doppler-broadened line is given by 



<l>y(z/ - Z/q) 



1 Az/p 
7r3/2 AZ/I5 



1 



exp 

5 



z/ — z/ 



/\2 



where Az/p and Az/^) meas ure the pressure and Doppler broadening widths, respective ly, and 
$y is normalized to unity ( Goody fc Yung 1989 : Thomas fc Stamnes 2002; Lioull2002l ) . Pres- 
sure broadening generally depends on pressure, temperature and the radiative constituent 
under consideration but the "classical" scaling 



Az/, = 0.02 - 0.05 



P 



T 



-1/2 



cm 



-1 



(9) 



1 bar / V 1500 K 

with a —1/2 power l aw dependence on temperature, should be sufficient for our order-of 
magn itude estimates ( Goody fc Yung 1989 : Thomas fc Stamnes 2002 : LioiJl2002 : Sharp fc Burrows 
20071 ). Doppler broadening depends on temperature, wavelength (Aq = c/z/q) and the mass 
mmoh of the radiative constituent under consideration. 



Az/^ ^ ^ ^ ^ l^^cs ^ 14 i ^ ^ ' cm-' 

^ C V mmol cV7"^mol ' ' V"^mol/ Vl^OOKy \1 filR J 

(10) 

where is the mass of the H2 molecule (the dominant atmospheric constituent), 7 is 
the gas adiabatic index and Cs is the corresponding adiabatic sound speed. This scaling 
illustrates how the Doppler width is reduced for radiative constituents which are typically 
more massive than molecular hydrogen. Note that the linear scaling of Aud with the central 



''The photosphere can be defined as the height in the atmosphere at which photon escape to space becomes 
possible, i.e., where the monochromatic optical thickness approaches unity. 
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wavelength Aq implies significant Doppler width variations across the relevant optical-IR 
spectral range. 

Substantial variations in temperature on constant pressure levels ar e found, from day- to 
night - side, i n current atmospheric circu l ation models for hot Jupit e rs (ICooper fc Showman 



20061: iLangton fc Laughlin 



200i 



2007 



Dobbs-Dixon fc Lin 2008: Showman et al. 



2008a 



^ Rauscher &: Menou 2009 ). Here, for simplicity, we choose 



20091 : iMenou fc R,auschei 
T = 1500 K and 1000 K as represent ative temperature values at the 1 bar and 10~^ bar 
levels, respectively (see, e.g.. Fig. 3 of iRauscher &: Menoul l2009l ). The ratio of Doppler to 
pressure broadening widths is given by 



3-7 



mH2_ 



1/2 



T 



1500 K 



P 



1 bar 



Ar 



1 yum 



:iii 



For relevant molecules, typical values of the mass ratio factor, (mHa/'^moi)^''^, are ~ 
1/2.8 (CH4), 1/3 (H2O), 1/4.7 (CO2) and 1/5.6 (TiO). At 1 bar, for T = 1500 K, the 
typical broadening ratio is thus Auu/Aup ~ 0.5 — 2 at 1 /im and 0.05 — 0.2 at 10 /xm. At 
10~^ bar and T = 1000 K, the ratio becomes Aud/Aup ^ 30 — 130 at 1 /xm and 3 — 13 at 
10 /xm. Doppler broadening is thus significant at the 1 bar level, especially in the near-IR 
and the optical, and it becomes increasingly dominant across the entire optical-IR spectral 
range higher up in the atmosphere. 

The sizable contribution of thermal Doppler broadening to the width of radiative lines 
at photospheric levels in hot Jupiter atmospheres, together with the scaling Aud < CgX uq/c 
from Eq. (fTOj) . suggests that bulk Doppler shifts from atmospheric motions near the sound 
speed could modify the shapes of radiative lines significantly. By contrast, much deeper 
in the atmosphere, where Auo/Aup <^ 1, Doppler shifts from bulk atmospheric motions 
near the sound speed would only amount to small shifts over comparatively wide, pressure- 
broadened lines. Furthermore, wind speeds themselves may be reduced at these deeper levels 
(see, e.g., FigJI]). 

These qualitative arguments are not very informative about the possible consequences 
of bulk Doppler shifts on radiation transport in a dynamic atmosphere. In particular, since 
the Doppler cores of radiative lines are often very optically thick ( "saturated" ) in hot Jupiter 
and other planetary atmospheres, they do not necessarily contribute much to the overall at- 
mospheric energy budget, at least in the static case0 To help us evaluate more quantitatively 



^Although this is not always explicitly stated, line-by-line radiative transfer models for hot Jupiter 
atmospheres currently available in the literature, with typical spectral resolutions ^ 1 cm~^ (e.g., 



the magnitude of bulk Doppler shift effects, we now turn to models of line transmissivities 
and equivalent widths in dynamic atmospheres. 



4. Transmissivity and Equivalent Width Models 

4.1. Transmissivities 

As summarized in ^ monochromatic transmissivities encapsulate the absorption and 
emission properties of the atmospheric medium (eqs. [2H3]). Here, we isolate the effects of 
velocity gradients along an arbitrary optical path by modeling the monochromatic transmis- 
sivity in a dynamic atmosphere as 



% = exp 



^0, 



^ I'proj 



X s') 



ds' 



:i2i 



where s' is the length along the unit optical pathj^ the line integrated opacity fctot is assumed 
to be constant along the path and is the dimensionless Voigt line profile defined in Eq. (|H]). 
The above expression for the bulk Doppler shift term in the line function assumes a constant 
velocity gradient along the path, i.e. a velocity offset that increases linearly with s', from 
at s' = to the maximum value Vproj at s' = 1. This simple model isolates the effects 
of bulk Doppler shifts by assuming that the path is otherwise homogeneous. In a more 
realistic atmospheric model, optical paths would be inhomogeneous, with line strengths, 
fctot, generally varying with temperature and line shapes, generally varying with both 
temperature and pressure along the specific path (iGoody &: Yunej|l989l : iThomas &: Stamnes 



20021 : lLioull2002[ ) 



pr 

H 



Sharp &: Burrows 2007 : Seager fc Sasselov IQQSt Seager et al. 2005h . do not necessarily resolve the narrow 



Doppler cores of radiative lines, with typical widths Az^d ^ 1 cm ^ according to Eq. (fTOl) . 

^In Eq. , fctot has unit of inverse length. Without loss of generality, we only consider optical paths of 
unit length in our simplified model. As a result, fctot fully characterizes the optical thickness of the modeled 
path. To differentiate the more general line integrated opacity /ctot appearing in Eq ([?]), which has units of 
cross-section per unit mass, from the simpler formulation adopted above, we write it as fctot in the simplified 
model. 

^In principle, the combination of Eqs. ([3]), ^ and ([S]) requires one to account for the profile variation 
with pressure and temperature separately for each individual line along the optical path of interest. A very 
common simplifying assumption in the treatment of atmospheric radiation, known as the Curtis-Godson 
approximation, is to treat the optical path as if it were homogeneous, i.e. with constant pressure and 
temperature, and use adequately path-averaged values for the line strength, shape and the absorber amount 
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Figure |2] shows representative profiles of monocliromatic transmissivities for a single line 
centered at u = uq in the model described by Eq. flT2|) . Fixed values for the line integrated 
opacity, fctot = 3, and for the ratio of Doppler to pressure broadening widths, Ai/£)/Az/p = 10, 
were adopted. The deepest transmissivity curve corresponds to the static reference model, 
with Vjproj = 0. In decreasing order of profile depth, from left to right, the other curves 
correspond to models with total velocity offsets Vproj = 1,3,5, 7 and 9 Auc x c/z/q, i.e., in 
units of the thermal Doppler width. 

From the point of view of absorption, transmissivities 71, < 1 correspond to the minimal 
absorption of the optically thin regime, while 71- <^ 1 correspond to the strong absorption 
limit of the optically thick regime. As is clear from Figure [2l monochromatic transmissivities 
can be considerably affected by bulk Doppler shifts approaching or exceeding the thermal 
Doppler width, Aud, of the line under consideration. A switch from partially optically-thick 
to fully optically thin occurs in Figure |2l While only the case with /ctot = 3 is shown, the 
shift and the flattening of line transmissivities seen here is qualitatively representative of 
what is seen for more general cases, with values of fctot ranging from ^ 1 to ^ 1. We note 
that the fiat-top nature of line transmissivities shown in Figure [2] when velocit y gradients are 



large, i.e., Vproj ^ x c/z/q, is consistent with the flattening discussed by ICastorl (119701 ) 
using the Sobolev approximation (see his Fig. 4). 

It is worth emphasizing that the magnitude of the velocity gradients considered in 
Figure [2] are relevant to hot Jupiter atmospheres. Writing the bulk Doppler shift term in 
Equation ([T]) as Az^buik = 'n-v{r) x z/q/c, we can express its overall magnitude in our simple 
transmissivity model in terms of the thermal Doppler width, Aud, as 



Az/buik 
Aud 



Tmmoi 
2m 



proj 



H2 



(13) 



Since, as seen earlier, the square root factor can reach values up to 5-6 for relevant 
molecules and, as suggested by Figure [H differentials of projected velocities can reach up 
to Vproj ~ 1^2 Cs along some atmospheric optical paths, the entire range of velocity offset 
values considered in Figure [2] is probably of interest for hot Jupiter atmospheres. 



(e.g., iRodgers fc Walshawl Il96d : ICoodv fc llQSOl : iThomas fc Stamnesl bood : IlIouI l2002h . While the 

accuracy of the Curtis-Godson approximation has been extensively tested in the context of static atmospheric 
radiation transport, its possible breakdown when optical paths acquire anisotropic properties in the presence 
of bulk Doppler shifts could lead to subtle complications in the treatment of radiation transport in a dynamic 
atmosphere. This specific aspect of the problem is not explicitly addressed by our simple transmissivity 
models. 
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By themselves, offsets and distortions of line transmissivities of the type shown in Fig- 
ure [2] are not observationally too significant because they occur over spectral intervals typ- 
ically < 1 cm~^, which is well below the effective resolving power achieved in optical-IR 
spectroscopic studies of hot Jupiter atmospheres. To the extent that such bulk Doppler 
distortions can alter the radiative energy balance in the dynamic atmosphere, however, they 
could in principle modify the radiative structure of hot Jupiter atmospheres. To evaluate 
the magnitude of this effect, we turn to a discussion of line equivalent widths. 



4.2. Equivalent Widths 



Let us first justify our use of equ ivalent widths by borr o wing from the existing literature 



2002: Liou 



on " narrow-band" spectral models (iGoody fc Yung] Il989l : iThomas fc Stamnes 
2OO2I ) . Ignoring the boundary term for simplicity, the integral term for the intensity given 
in Eq. ([2]) is essentially of the form 



B^dTu- 



(14) 



which is simply reformulating it as the cumulative Planck emission and its absorption along 
a specific optical path, weighed by the transmissivities of the various contributing layers, 
from closely adjacent ones (with transmissivity 7^ ~ 1) to more distant ones (with 71, — 0, 
in the optically thick limit). In narrow-band models, the frequency- integrated intensity is 
written 



/•oo i*\ n\ 

Jo Jo -'0 



(15) 



where Az/j is the frequency span of the zth narrow band, the Planck function value Bi is 
considered to be a constant over each narrow band and the band-averaged transmissivities 
are defined by 



Tvdu. 



(16) 



The narrow-band formalism is well justified as long as the number of narrow bands is 
sufficiently large for the Planck function to be we l l approxim ated by a constant in each band 
(iGoody fc Yung||l989l : iThomas fc Stamned l2002l : lLioull2002l ). Since atmospheric fluxes are 
obtained by angular integration of the intensity, equation f lTSj) relates the radiative energy 
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balance of the atmosphere to multiple Planck- weighted integrals of band-averaged transmis- 
sivities. 

The narrow-band formalism is traditionally used with bands that are still wide enough 
to encompass a large number of individual radiative lines. By contrast, we use it here in its 
simplest formulation, with very narrow bands including only a single s pectral line, assuming 
that each such line is well isolated from all the other lines (see, e.g., iGoody &: Yung)ll989 



chapter 4). In that limit, the equivalent width (EW) of a specific line can be written 



EW= I {l-%)dv = Avi{l-Ti), (17) 

where Az/j now represents the typical spacing between neighboring isolated lines. It should be 
noted that the EW is independent of the value of Ai^j adopted as long as the spectral interval 
is large enough to encompass all of the transmissivity values contributing meaningfully to 
the above integral (i.e., % < I). Thus, when line overlap can be omitted, Eqs. f fT5] - [T7|) 
clarify the direct relation that exists between the radiative energy balance of an atmosphere 
(oc /), Planck- weighted integrals of band-averaged transmissivities (oc dTl) and the equivalent 
widths of all important radiative lines {dTi oc —dEW). This is our main justification for 
using deviations in line equivalent widths as a quantitative measure of the effects of bulk 
Doppler shifts on the overall radiative balance of a dynamic atmosphere. 

We perform equivalent width calculations by integrating transmissivity profiles like the 
ones shown in Figured], following equation f lT7|) . For the calculations presented here, we 
verified that our EW results are independent of the spectral interval chosen for integration, 
as long as it spans more than 100-1000 thermal Doppler widths (Az/^i) on each side of the line 
central frequency, uq. We evaluate the effects of bulk Doppler shifts simply by comparing the 
equivalent widths of lines in a static atmosphere, EWq = EW{Vpj-oj = 0), to those obtained 
in the presence of finite bulk Doppler shifts (Vproj 7^ 0). 

In all our calculations, using the simple transmissivity model described by Eq. (fT2|) . 
we find that EWs in dynamic atmospheres are systematically larger than the corresponding 
value in a static atmosphere, EWq. We quantify this trend with the fractional excess, 
EW/EWq — 1, expressed in percents. For example, for the various shifted transmissivity 
profiles shown in Figure [2l we find EW excesses of 2.4, 14, 24, 31 and 35% over the static 
EWq value for the profiles with total velocity offsets Vp^oj = 1,3,5,7 and 9 Au^, x c/uq, 
respectively. 

Figure [3] further clarifies the variation of the EW excess with the magnitude of the 
total velocity offset, Vproj, in units of Auo x c/z/q, in a model with fixed ratio of broadening 
widths, Auu/Aup = 10, for the Voigt line function. The various curves show models with 
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line integrated absorption coefficient fctot = 0.3 (dashed line), 3 (dotted; same as in Fig. |2]), 
30 (solid), 300 (dash-dotted) and 3000 (triple-dot dashed). While the EW excesses increase 
with the value of Vproj, the exact dependence is not entirely trivial. For k^ot = 3, 30 and 300, 
EW excesses from several tens of percents to more than a hundred percent are possible, for 
the same range of Vproj values as deemed relevant for hot Jupiter atmospheres earlier. 

Figure H] shows the systematic variation of the EW excess with the magnitude of the 
total line absorption coefficient, ktot, in a model with a value of the total velocity offset fixed 
at Vpj-oj = 10 Aud X c/z/q. The various curves show models with Voigt function broadening 
ratios Auu/Aup = 1 (dashed line), 10 (dotted; same as in Figs. |2]and[3]), 100 (solid) and 
1000 (dash-dotted). The EW excess ffist increases with the value of fctot, peaks around 
ktot ~ 20-40 and then drops at larger fctot values. Peak EW excesses are thus reached for 
rather strong optically thick conditions. Note that EW excesses ~ 150-200% are possible 
for ktot ~ 30, Auo/Aup > 100 and Vproj ~ 10 Au^ x c/uq. EW excesses easily reach several 
tens of percents at fctot ~ 30 in all of the models considered here, for the rather large value 
of V^roj = 10 Aud X c/z/q adopted (see again Fig. [3] for the dependence with Vproj)- 

We can understand various trends in the behavior of the EW excesses with V^roj, ^tot 
and Aud/ Aup in our simple models as follows. As already mentioned earlier, we have found 
that the shifts and distortions of monochromatic transmissivity proffies shown in Fig. [2] are 
qualitatively representative of the general behavior seen in all our models, whether the line is 
in the optically thin or optically thick regime. For a line well into the optically thick regime 
{1 — — 1), the growth in EW with l^roj can be simply understood as a stretching of 
the optically thick portion of the transmissivity curve, which contributes maximally to the 
growth in EW according to the integral in equation (fT7|l . This is the main trend observed 
for moderate to high values of fctot in Figure |3l 

In the optically thin regime, the increase in EW excess with ktot can be understood 
as resulting from the presence of larger radiative intensities outside the line center (where 
maximum absorption occurs), so that Doppler-shifted absorption around the line center 
contributes increasingly to the EW excess. As the value of fctot is increased well into the 
optically thick regime, however, the much shallower pressure-broadened wings of the Voigt 
line profile start making a significant contribution to the EW integral, over an increasingly 
wider range of frequencies. This reduces the effective contribution of the optically-thick, 
shifted Doppler core to the total line EW, which explains the decline in EW excess at large 
ktot values in Figure HI the more so in models with larger values of the pressure-broadening 
parameter, Aup. 
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5. Discussion and Conclusion 

Our main result has been to establish the possibility that bulk Doppler shifts from fast, 
sonic or transonic wind speeds in hot Jupiter atmospheres could affect radiation transport 
and thus in principle influence the overall radiative structure of these atmospheres. We have 
quantified the magnitude of these effects with simplified transmissivity models in a dynamic 
atmosphere, by assuming that linear bulk velocity gradients are present along otherwise 
homogeneous optical paths and by adopting Voigt line profile parameters appropriate for 
the radiatively forced regions of hot Jupiter atmospheres. We have found that bulk Doppler 
shifts systematically increase the equivalent widths of isolated radiative lines, relative to 
the equivalent widths in a static atmosphere, with excesses easily reaching several tens of 
percents for relevant model parameters. 

Flux errors below the few percent level may well be acceptable in atmospheric models 
currently used to interpret the increasingly rich set of observational data on hot Jupiter 
atmospheres since other important atmospheric unknowns also contribute to modeling errors 
(e.g., the exact compositional profiles of radiatively active constituents). Flux errors at the 
level of several tens of percents or more may be too large to be ignored, however, especially if 
they originate from systematic excesses in the equivalent widths of important radiative lines. 
Given the link between line equivalent widths and radiative intensities discussed in §4.2l in the 
context of the narrow-band formalism, we therefore suggest that equivalent width excesses 
of the magnitude found in several of our idealized models could have a significant impact on 
radiative fluxes and the overall radiative energy balance of hot Jupiter atmospheres. 

However, we must also caution that, beyond these preliminary arguments, it is not pos- 
sible to settle this question without a considerably more detailed calculation than the one 
presented here. Indeed, it is the combined effect of bulk Doppler shifts on a vast number of 
radiative lines, all with varying degrees of optical thickness, weighted by the various Planck 
functions of the atmospheric layers under consideration (thus depending on the atmospheric 
temperature profile itself), along a variety of inhomogeneous optical paths, which ultimately 
determines by how much the overall radiative balance of the dynamic atmosphere is af- 
fected. By contrast, our simple transmissivity and equivalent width models were focused 
on single isolated lines with purely linear velocity gradients along arbitrary, homogeneous 
optical paths. Furthermore, contributions from continuum opacity sources in the atmo- 
sphere, for instance in cloudy regions, could reduce the impact of bulk Doppler shifts on 
the atmospheric radiation transport, since small shifts mostly affect narrow radiative lines. 
Various simplifying assumptions made in our models would thus be invalidated in a more 
realistic radiative transport calculation. Interestingly, the full radiation transport problem 
in a dynamic atmosphere may be amenable to practical numerical solutions in the future 
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e.g., knop et all 120091 ) 



It is also worth emphasizing that our discussion has largely focused on the long- wavelength 
(thermal) component of the atmospheric radiation spectrum, for which the diffuse approxi- 
mation and the use of an isotropic (Planck) source function are justified. At first, one may 
be tempted to neglect the effects of bulk Doppler shifts from horizontal winds in the treat- 
ment of short-wavelength atmospheric radiation since the corresponding beam-like insolation 
is vertical in first approximation. This could be misleading, however. When scattering is 
important, as is usually the case, it generates optical paths slanted enough that they could 
become much more susceptible to the bulk Doppler shifts caused by horizontal winds. Fur- 
thermore, there is a significant fraction of the atmosphere that is located far enough away 
from the substellar point to be subject to partially slanted irradiation. In the case of close-in 
planets like hot Jupiters, the finite size of the stellar disk could also contribute to slanted 
irradiation and thus an increased sensitivity to Doppler shifts from horizontal winds. Let 
us emphasize that radiation transport in a dynamic atmosphere with an anisotropic source 
function, as would result from multiple scattering of the beam-like insolation in the short- 
wavelength portion of the atmospheric radiation spectrum, would lead to a considerably 
more complex formulation of the radiative transfer problem than discussed here for the case 
of diffuse thermal radiation (e.g., Eq. [2]). For example, one wonders whether the anisotropic 
transport resulting from atmospheric bulk motions could result in a stronger polarization 
signal than has b een estimated on the basis of static radiative transfer calculations (e.g.. 



Seager et al.ll2000f ) 



Interestingly, the various short-wavelength radiation transport effects mentioned above 
should be particularly pronounced in the context of transmission spectroscopic measure- 
ments, which specifically probe nearly horizontal optical paths at the atmospheric planetary 
limb during transits. Essentially all transit spectrosc opic di a gnost ics could thus be affected 



by bulk Doppler shifts from fast atmospheric winds. iBrownl (120011 ) has presented a detailed 
discussion of this problem, including consequences of the vertical gradients of horizontal wind 
velocity. The possibly more important effects of bulk Doppler shifts along a specific optical 
path were omitted from these calculations, however, even though the author did comment 
on the expectation of increased line equivalent widths. It may thus prove important to 
carefully reassess various interpretations about the chemical composition of hot Jupiter at- 
mospheres based on transit spectroscopic measurements with models which properly account 
for radiation transport in dynamic, rather than static, atmospheres. 

Finally, let us conclude by mentioning a few possible extensions of this work to other 
classes of planetary atmospheres. It would seem that the standard static assumption for 
the treatment of radiation transport in Solar System planetary atmospheres is well justified. 
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Indeed, for these much cooler atmospheres, at a representative pressure level of 1 bar, radia- 
tive lines are significantly more pressure-broadened in Solar System atmospheres than in hot 
Jupiter atmospheres (see scalings in §3.2p . In that limit, and for wind speeds well below the 
sound speed, our models do indicate very small excesses in the equivalent width of radiative 
linesj^ By contrast, the typically much hotter planets discovered by astronomers in recent 
years, which are often subject to unusually strong radiative forcing conditions, would seem to 
be more natural sites for the application of the radiation-hydrodynamical principles empha- 
sized in the present work. Besides hot Jupiters, these principles could also find applications 
in the class of eccentric giant planets which experience t r ansient atmospheric flash-heating 



during periastron passage (e.g.. lLangton &: Laughlinll2008l : iLaughlin et al.ll2009[ ) and perhaps 



the emerging class of hot super-Earths, if the atmospheres of these planets are able to sus- 
tain wind velocities approaching or exceeding the sound speed in the presence of significant 
ground drag. 
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Fig. 1. — Representative profiles of wind velocities, in units of the local sound speed c^, 
projected along various opti cal paths at a fixed z enith angle 6*0 = 53 degrees in the atmo- 
spheric circulation model of iRauscher &: Menoul (120091 ). The various panels show profiles 
at the model substellar point (a), antistellar point (b), west equatorial terminator (c) and 
north pole (d). In each panel, the various curves show projected velocity profiles for optical 
paths oriented to the east (solid line), north-east (dotted), north (dashed) and north-west 
(dash-dotted). Velocity differentials which are a sizable fraction of, and in some cases exceed, 
Cs are typical. 
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Fig. 2. — Monochromatic transmissivity curves for a line centered ai u — i^q with a Voigt 
profile shape and a fixed ratio of Doppler- to pressure-broadening widths, Az/£)/Az/p = 10. 
From left to right, transmissivity profiles are shown for linear bulk velocity gradients with 
total offsets Vproj/^i^n = (static), 1,3,5,7 and 9 xc/uq over the full optical path. The 
specific model shown here is marginally optically-thick, with a line-integrated absorption 
coefficient /ctot — 3 over the full optical path, but these transmissivity profiles are represen- 
tative of more general cases. Monochromatic transmissivities are significantly affected by 
bulk Doppler shifts in excess of a few thermal Doppler widths. The increasingly blue-shifted 
transmissivity profiles shown have equivalent widths exceeding that of the static profile by 
2.4, 14, 24, 31 and 35%, from left to right. 
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Fig. 3. — Equivalent width excess (in %) of a radiative line in a dynamic atmosphere with a 
linear bulk velocity gradient, as a function of Vpj-oj, the total velocity offset over the optical 
path, in units of the line's thermal Doppler width {Auu x c/uq). In this specific model, the 
line has a Voigt profile shape with a fixed ratio of Doppler- to pressure-broadening widths, 
Aud/Aup = 10. The various curves show results for fines with different total absorption 
coefficients: ^tot — 0.3 (dashed line), 3 (dotted), 30 (sofid), 300 (dash-dotted) and 3000 
(triple-dot dashed). 
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Fig. 4. — Equivalent width excess (in %) of a radiative line in a dynamic atmosphere with 
a linear bulk velocity gradient, as a function of the total line absorption coefficient, /ctot- In 
this specific model, a relatively large value of the total velocity offset over the full optical 
path, V^roj = 10 Ai/£, X c/^o, is adopted. The line has a Voigt profile shape and the various 
curves shown correspond to different ratios of the Doppler- to pressure-broadening widths: 
Ai/^/Ai/p = 1 (dashed hne), 10 (dotted), 100 (sohd) and 1000 (dash-dotted). 



